Exciton supersolidity in hybrid Bose-Fermi systems 
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We investigate the ground states of a Bose-Einstein condensate of indirect excitons coupled to an 
electron gas. We show that in a properly designed system, the crossing of a roton minimum into the 
negative energy domain can result in the appearance of the supersolid phase, characterized by peri- 
odicity in both real and reciprocal space. Accounting for the spin-dependent exchange interaction of 
excitons we obtain ferromagnetic supersolid domains. The Fourier spectra of excitations of weakly 
perturbed supersolids show pronounced diffraction maxima which may be detected experimentally. 
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A supersolid phase is characterised by superfludity and 
long range order in both real and reciprocal spaces [1|. 
One of the suggested indications of this phase transition 
is a "rotonlike instability" observed if the energy of the 
roton minimum of a superfluid is pushed below the energy 
of the ground state (characterised by zero wave- vector, 
q) 0. Strong indications of the phase tranistion to su- 
persolidity have been found in 4 He at low temperatures 
and high pressure [3]. Until now, experimental stud- 
ies of supersolidity have been limited to 4 He and 3 He, 
while, potentially, this phase may occur in any system 
of interacting bosons [4] . Here we discuss the possibility 
of formation of exciton supersolids in specially designed 
semiconductor structures. 

Recently, Bose-Einstein condensation (BEC) and su- 
perfluidity of crystal quasiparticle excitations - excitons 
or exciton polaritons - have been reported in various 
semiconductor systems [1-0 • The exciton-exciton inter- 
action is usually almost independent of the wave vector 
Q, which is probably why in most experimentally stud- 
ied cases the roton minimum is not formed. Neverthe- 
less, formation of a roton minimum and even of a roton 
instability has been predicted for a hybrid system where 
a two-dimensional (2D) exciton condensate is placed in 
close vicinity to a 2D electron gas [Toj| . Moreover, the 
recent theoretical paper by Parish et al [ll| evokes the 
possibility of the supersolid phase for excitons in asymet- 
rically populated bilayers. Here we show that an exciton 
condensate in a hybrid system where excitons and an 
electron gas interact but do not overlap can undergo a 
transition to the supersolid phase. 

We shall consider a system of coupled semiconductor 
quantum wells (CQW) which contain the BEC of spa- 
tially indirect excitons which has been widely studied in 
the past [12]. We introduce an important new element 
to this system: a 2D gas of free electrons located in a 
thin metallic layer grown on the top of semiconductor 
structure. The metal is separated from CQW by a wide 
band gap semiconductor layer, so that there is no over- 
lap between free carriers and excitons, while they are 



close enough to each other to allow for efficient Coulomb 
interaction. The system is biased in such a way that 
the exciton dipole moments are oriented towards metal. 
This configuration prevents attraction between excitons 
and electrons and formation of trions. We phenomeno- 
logically describe spin-dependent interactions of dipole 
excitons by two constants, as specified below, and con- 
centrate on the new physics brought by exciton-electron 
interactions. The system we study has been proposed for 
realization of charged electron-hole complexes [13] and 
of superconductivity mediated by the exciton conden- 
sate [14[ . Compared to the unequally populated electron- 
hole bilayers considered by Parish et al pjj this system 
has an evident advantage of protecting excitons from the 
phase space filling effect, which would allow operation 
with much higher exciton and 2D electron gas densities 
than in any system where excitons and free carriers spa- 
tially overlap. In the present work we concentrate on the 
phase trasitions of the exciton condensate induced by the 
presence of an electron gas in the neighboring metallic 
layer. Referring to the most popular GaAs based system, 
we account for the exchange interaction between excitons 
having +1 and —1 spin projections to the structure axis. 
This interaction, which represents the superexchange via 
intermediate dark excitons with spin projections +2 and 
—2, has been recently addressed theoretically TH, 16| . 

An exciton condensate is described by a fac- 
torized Hartree-Fock ansatz ^±(ri, . . . , rjv, t) = 
Yli 1 P±( T iit) [17], which, in the limit of large number 
of particles N, leads to the following mean field energy 
with a nonlocal interaction potential, where ± denotes 
the +1 and —1 exciton spin states, 
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FIG. 1. The excitation spectrum (top) and effective exciton- 
exciton interaction potential (bottom) for three values of 
the exciton-electron layer separation L. We use the exci- 
ton mass m ex = 0.21ra e , and the average exciton density 
ntot = 10 11 cm -2 . The exciton dipole separation, / = 12nm, 
and the correction factor, £ = 0.05. 



N + \^ + (r,t)\ 2 + JV_|^_(r,t)| 2 , and V e f_ ex (r) is the spin- 
isotropic exciton-exciton interaction potential. The cor- 
responding time-dependent Gross-Pitaevskii (GP) equa- 
tions for the two condensates read 
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tion depending on frequency and wave- vector reads 
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with e denoting the dielectric constant of the media and A 
the sample area. The matrix element of exciton-exciton 
interaction V22 is dominated by dipole-dipole repulsion 
[H EH and can be expressed as 
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where d is the exciton dipole moment in normal to the 
QW plane direction, £ accounts for the reduction of 
dipole-dipole repulsion due to electron-hole pair corre- 
lations in the exciton condensate [19]. The expression 
for the matrix element of electron-exciton interaction 
reads (the details of the calculation can be found in e.g. 
RefH)j): 
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where /3 e ^ = m e? ^/(m e -f ra^), with me,/! being the effec- 
tive masses of electrons and holes. 

The poles of the effective potential determine the dis- 
persions of the collective modes of the exciton system 
given by the equation 

(M 2 = (E ex (q)) 2 + 2N V 22 (q) E ex (q), (7) 

where 

V 22 (q)=V 2 2 + V 1 %(q)n 1 (q,co)/[l-V 11 (q)U 1 (q,co)}. (8) 



Derivation of the effective exciton-exciton interaction In general, the exciton-exciton interaction mediated by 
potential in this system has been recently published [Io|. the electron gas is a time-dependent, retarded inter ac- 
The matrix element of effective exciton-exciton interac- tion. In order to obtain a non-retarded GP equation (j2]), 
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FIG. 2. Examples of density profiles of the exciton condensate 
in real space (left column) and reciprocal space (right column, 
density in logarithmic scale). The figures (a), (b) show the 
homogeneous ground state for separation L = 30 nm, and 
the figures (c), (d) depict the supersolid ground state for 
L = 18 nm, characterized by periodicity in both real and mo- 
mentum spaces. The orientation of the supersolid lattice is 
chosen spontaneously, (e), (f) A metastable state with lattice 
imperfections, corresponding to a local energy minimum. In 
momentum space, such a state is characterized by concentric 
rings. 



we average the potential (|3j) over a frequency range lim- 
ited by the exciton dissociation frequency Ud = Eb/Ji. 
Depending on the frequency- averaged interaction poten- 
tial V2jjl ex (q), the excitation spectrum of the condensate 
coupled to an electron gas can develop a roton minimum 
due to the attractive character of the electron mediated 
interactions [lOj, see Fig. [TJ This Figure shows the exci- 
ton dispersion curves calculated from Eq. ([7]) in compar- 
ison with the average interaction potential for different 
separations L between the QWs containing the 2D elec- 
tron gas and the exciton condensate. Indeed, for small 
separation between excitons and electrons the potential 
becomes attractive in an intermediate range of q and a 
roton instability develops. 

The existence of a roton instability has been connected 
with the possible existence of a supersolid phase [l], Q , 



in which superfluid and crystalline orders coexist. We 
investigate this possibility by numerical minimization of 
the Hamiltonian (pQ) using the imaginary time propaga- 
tion method. To start with we assume a = and in- 
vestigate a single-component exciton condensate. We de- 
termine the ground states for various physical parame- 
ters in the case of a two dimensional box with periodic 
boundary conditions. We take realistic values of param- 
eters m ex = 0.21ra e , m e being the free electron mass, the 
average exciton density n to t = 10 n cm -2 , with an exci- 
ton dipole separation I = 12 nm, and different effective 
interaction potentials corresponding to various exciton- 
electron layer separation L. 

We distinguish three regimes in which the minimiza- 
tion procedure leads to qualitatively different results. In 
the case of a positive roton gap, when the excitation spec- 
trum does not cross into the negative energy domain, the 
ground state is a homogeneous superfluid, as depicted 
in Figs. [2fa) and (b) for L = 30 nm. When decreasing 
the interlayer distance to L = 18 nm, a dramatic change 
occurs, and a state characterized by a periodic triangu- 
lar lattice becomes the ground state of the system, see 
Figs.[2]^c) and (d). The orientation of the lattice is chosen 
randomly in a process of spontaneous symmetry break- 
ing, depending on the initial fluctuations in the numerical 
minimization procedure. This symmetry breaking is ad- 
ditional to the symmetry breaking associated with the 
randomly chosen phase of the condensate wave function. 
In the reciprocal space, the appearance of crystalline or- 
der is manifested by the appearance of side peaks corre- 
sponding to the lattice wave vectors. We note that the 
lattice constant of the pattern is determined by the wave 
vector value at the roton minima. Indeed, the side peaks 
closest to k — in Fig. Efd) match with the position of 
the roton minimum in Fig. [TJ 

When decreasing the interlayer distance even further 
to L = 14 nm, the minimization algorithm leads to the 
collapse of the wave function ip± into a single point of 
the numerical mesh. This indicates that the ground state 
cannot be determined due to the breakdown of the mean- 
field model. To obtain a correct ground state, one has to 
include higher order effects that would limit the growth 
of the local condensate density. We note that collapse is 
generally encountered for a nonlinear Schrodinger equa- 
tion in the 2D (critical) case with attractive nonlinear- 
ity [2i|. 

In the supersolid regime (L = 18 nm), depending on 
the initial noise, the minimization procedure can produce 
yet another type of solution, which resembles a crystal 
structure with dislocations, see Fig. Efe). While these 
solutions do not correspond to the ground state of the 
system, they are local minima of the energy and consti- 
tute metastable states, and as such are likely to appear 
in experimental realizations. In the momentum space, 
these disordered supersolid patterns are characterized by 
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FIG. 3. Numerically determined spectra of a weakly per- 
turbed (a) noninteracting exciton condensate, (b) condensate 
interacting through a local (momentum-independent) poten- 
tial, and (c) supersolid condensate in the presence of coupling 
with electron gas, interacting through an effective potential 
Kx^ex (parameters as in Figs.[2fc),(d)). While in the cases (a) 
and (b) we obtain the parabolic and linear Bogoliubov spec- 
tra as predicted by the theory, in the case (c) the spectrum is 
distorted as a result of interference of branches starting from 
different /c-components of the ground state in the momen- 
tum space, corresponding to the vertical lines. The density is 
shown in a logarithmic intensity scale. 



concentric rings around the k = point, see Fig. EJf). 

We also investigate the influence of the supersolid tran- 
sition on a weakly excited system. In Fig. [3] we show the 
spectra of perturbed ground states allowed to evolve in 
time according to Eq. (j2j). We obtain the spectra by 
Fourier transforming the wave function ^ ex in space and 
time. The three figures correspond to the cases of a non- 
interacting exciton condensate, a condensate with a q- 
independent interaction potential, equivalent to a contact 
(local) repulsive interaction in the absence of an electron 
gas, and an exciton supersolid interacting through a q- 
dependent interaction potential V^L ex . While the first 
two figures show the corresponding free-particle parabolic 
spectrum and the Bogoliubov spectrum of excitations, 
respectively, in the case of the supersolid the spectrum 
becomes blurred as a result of interference of several 
branches starting from different points in the momentum 
space. These points correspond to the /c-space compo- 
nents of the ground state, visible as pronounced vertical 
lines. The excitation branches appear to be linear close 
to their points of origin, which is an indication of a non- 
vanishing superfluid fraction. 

Interesting new physics can be expected in two- 
component exciton condensates, where a ^ 0. Fig. |4] 
shows an example of the spatial distribution of excitons 
in the supersolid phase in this case. One can see that 




FIG. 4. Numerically calculated density plot of the exciton 
condensate in the supersolid phase, showing the ferromagnetic 
domain structure. Here n+ — n_ — 10 n cm -2 , L — 18 nm 
and a is chosen in such a way that for q — the interaction 
of excitons with opposite spin is 10% stronger than the inter- 
action of excitons with parallel spins. The green and orange 
coloring corresponds to the +1 and —1 spin components. 



ferromagetic domains are formed: the triangular crystal 
lattice is maintained in each domain, while long-range 
order is broken by stochastic domain boundaries. 

Experimentally, formation of the exciton supersolid 
may be detected by angular and polarisation resolved 
photoluminescence, which would give access to the char- 
acteristic diffraction pattern (most likely, concentric 
rings) in the Fourier spectrum of excitations of the su- 
persolid. 
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